Spectral properties of the small polaron 
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We compute exactly both the spectral function of the electron and of the small polaron for 
the two site Holstein model. We find that for intermediary coupling, the small polaron is a better 
fundamental excitation of the system than the electron. However, the Lang-Firsov approximation 
fails to predict the right dispersion relation for the small polaron. 



Electronic properties of strongly interacting electron- 
phonon systems are now widely studied since it has been 
realized experimently that the formation of small po- 
larons in some of the new materials, such as the high 
Tc cuprates [|| or the perovskite manganates [||, is re- 
sponsible in part of some new and unusual physics. 

It is generally believed that the simplest model for the 
study of small polarons, is the so called Holstein molec- 
ular crystal model [||. In the weak coupling limit, elec- 
trons behave as good quasiparticles, while in the strong 
coupling limit the model is believed to exhibit well de- 
fined small polarons 0|. Perturbation theory towards 
intermediary coupling can then be used from the strong 
limit fixed point in term of small polaron quasiparticles 
1^,0. The validity of such a perturbation theory rely 
on the quality of the small polaron quasiparticle. If it 
is destroyed away from infinite coupling then we can no 
longer described the physics in term of small polarons 
and no longer use perturbation theory. 

In this paper, we compute exactly the spectral func- 
tions for the small polarons in two-site Holstein model for 
intermediary coupling. We find that they are good quasi- 
particles with a well defined relation dispersion w^. We 
can then compute the effective bandwidth of the small 
polarons D* = maa;(w^) — min{uj'l,) and show that it 
dramatically decreases with the strength of the electron- 
phonon coupling constant. These results clearly show 
that the physical picture given by Lang and Firsov is the 
right one. 

Our calculations are based on a two site cluster where 
exact diagonalization are easily performed. The hamil- 
tonian for the Holstein model is 



(1) 



The sum over S run on the nearest neighbours, t is the 
hopping integral for the tight binding approximation, luq 
is the optical frequency of the phonon and g is a di- 
mensionless coupling constant. The polaronic energy is 
Ep — g^ujf) and the usual elect ron-phonon coupling con- 
stant is A = Ep/zt where z is the number of nearest 
neighbours {z — 1 for the two site cluster). 



In the atomic limit, t = Q, the hamiltonian can be 
diagonalized exactly by the Lang-Firsov transformation 
The new hamiltonian is -ff = e^^He^ with S = 



bj)nj and rij 



Cj^a — UjCj^cr and bj = bj 



grij with 



The transformed hamiltonian is 



We thus obtain 



(2) 



- ^0 



J2b]b,-2EpY,nj, 



(3) 



The hopping term is now a complicated operator. In 
the Lang-Firsov approximation, one obtains an effective 
hamiltonian for the electrons by eliminating the phonon 
states, keeping only the vacuum without phonon. At zero 



temperature one gets Uj 



^ and thus an effective 



hopping integral t* = e^^ t for the small polaron. In the 
strong coupling limit, g — > oo, this approximation gives 
localised polarons, while for finite coupling it gives a pola- 
ronic band with the dispersion relation — —Ep-\-t*S,k, 
where is the Fourier transform of the kinetic term. 
For the two site cluster, one gets e'l = —Ep ± t* for 
k = 0,7T. If one considers systems with less than two 
electrons, then the interaction term responsible of the 
formation of bipolarons (the last term of eq.(^) vanishes 
and the small polarons energies ej^ are the right excita- 
tions of the system. The spectral function of a small po- 
laron, always in the Lang-Firsov approximation, is sim- 
ply ^(k, uj) = 2tt6{u! — £^). Inversely, the electrons are no 
longer good excitations. They carry a cloud of phonons, 
with an average number of phonon ng^, where n is the 
number of electron. The spectral function is given by 



A{k,Lu) = e-s 2ttS{u;-sI) 



1 

M 



(4) 



where M is the number of sites. Exact diagonalizations of 
small clusters |^ based on the bare hamiltonian of eq. (|l|) 
showed that the spectral function for the electron is too 
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much complicated to be interpreted (as can be seen in 
fig|^). For intermediary and strong coupling, there is no 
longer a quasiparticle peak and thus impossible to extract 
a dispersion relation. One possible test for the validity 
of the Lang-Firsov approximation should be to fit the 
numerical results for the spectral function of the electron 
with the formula of eq.(^). A more efficient method con- 
sist to compute directly the spectral function of the small 
polaron based on the transformed hamiltonian of eq.(|^). 
If the small polaron defined by the canonical transfor- 
mation is a good quasiparticle then the spectral function 
should contains only one main peak with a width smaller 
than the excitation energy. It is then easier to compare 
this energy to the one predicts by the Lang-Firsov ap- 
proximation. However the aim of this paper is simply to 
set up if the small polaron is a good quasiparticle for the 
system. 

In the transformed hamiltonian H, the operators Cj.o- 
and Uj commute. One can build separately the opera- 
tor cj+^^Cj^o- and then work in a subspace with a fixed 
number of polarons. The operator Uj is computed via its 
matrix elements given by 



n n 
p=0 p'=0 



{P+p') 
p\p'\ 



-,1/2 



[n — p)\{n' — p')\ 



(5) 



and the condition n — p ~ n' — p' . For the phonon basis 
states, a better convergence is obtained if we define some 
subspaces with constant total number of phonons . We 
diagonalize matrices with a maximum of 35 phonons (cor- 
responding to 666 basis states of phonons). Compari- 
son of the eigenvalues of both the bare and transformed 
hamiltonians H and H shows that the polaronic repre- 
sentation is better for the results presented here. In order 
to deal with the case wq — we choose the same param- 
eters as in ref. where the spectral function J{k,uj) 
for the electron has been computed, that is o^o = 1/1-1, 
g = 0.3-\/2, 0.8-\/2 and 1.3-\/2. We compute the spectral 
function at zero temperature in the ground state with 
zero small polaron defined by 



J{k,u) = 27r^ |(iV = l,m I 4 | TV = 0, 0) 



fe,T 



X 5{uj + E'^ 



N=0 



E. 



N=l 



(6) 



The number of phonon states has been chosen such that 
the results presented here are well converged. The cri- 
terium is that the transition amplitudes and energies con- 
verge towards a finite limit as we increase the number of 
phonon states. 
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FIG. 1. Spectral functions J{k = 0,tj) for the electron, 
for three different couplings and uq = 1/1.1, g = 0.3\/2 (a), 
5 = 0.8^ (b), g = 1.3^/2 (c). 

In fig.|], we show the spectral functions of the electron 
based on the hamiltonian of eq. (0) which correspond with 
the results of ref. |^ (we note that these are computed us- 
ing two different methods) and in fig]^ we show the spec- 
tral functions of the small polaron based on the hamil- 
tonian of eq.(^. The solid lines in figures correspond to 
Lorentzian of width O.lt while peaks correspond to the 
value of the transition amplitude, without the 27r fac- 
tor. We see that for weak enough coupling, both the 
small polaron and the electron are good quasiparticle 
with mainly one transition in the spectral function. As 
we increase the coupling we observe for the small polaron 
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that the spectral weight stays in one peak while it spreads 
in many peaks more or less spaced of the energy ujq for 
the electron. For g — 0.3-\/2 we find uJk=o — —1-108 
and iVk^T, ^ 0.977, for g = 0.8\/2, uJk=o ^ -1.795 and 
LUk=TT — —1.349, and for g — 1.3^/2, u!k=o — —3.314 and 
uJk=TT — —3.225. All the results are in unit of t. 
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FIG. 2. Spectral functions J(fc = 0,aj) for the small po- 
laron, for three different couplings and (jjq = 1/1-1, 5 = 0.3\/2 
(a), g = 0.8^2 (b), g = 1.3^ (c). 

If we compare the renormalized bandwidth D* ~ 



i^k=TT — i^k=o with the Lang-Firsov approximation D^p = 
2te~^ we find a discrepancy of 25% for the two strongest 
couplings. 

To conclude, the approach presented here is different 
from the standard approach of the so-called polaron prob- 
lem where one considers a system with only one electron 
coupled to the lattice and find the excitation energy of 
the whole Hamiltonian ||l^. Instead, we have shown 
that the small polaron defined through the canonical 
transformation of Lang-Firsov is a good quasiparticle for 
strong and intermediary electron-phonon coupling. This 
method applies directly to many body problems, such as 
the Hamiltonian of eq. (||) for the Holstein model and for 
arbitrary filling of the system. Perturbative expansion 
from e.g. the strong coupling Lang-Firsov approxima- 
tion fixed point MuM should be relevant. 
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